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Abstract 

The dynamics of Dvali-Gabadadze-Porrati Cosmology (DGP) braneworld with an anisotropic 
brane is studied. The Friedmann equations and their solutions are obtained for two branches of 
anisotropic DGP model. The late time behavior in DGP cosmology is examined in the presence of 
anisotropy which shows that universe enters a self-accelerating phase much later compared to the 
isotropic case. The acceleration conditions and slow-roll conditions for inflation are obtained. 
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I. INTRODUCTION 



Cosmological models inspired by higher dimensional theories have received a lot of at- 
tention in recent times. In these models our observed four-dimensional (4D) universe is a 
three dimensional hypersurface (brane) embedded in a higher dimensional space-time called 
bulk (for a review see jj). Two popular cosmological models emerging out of higher di- 
mensional theories are Randall-Sundrum (RS) [2] and the Dvali-Gabadadze-Porrati (DGP) 
model 3|. The generalisation of such models to the homogenious and isotropic Friedmann- 
Robert son- Walker (FRW) brane leads to modification of the Friedmann equation with a 
quadratic correction to energy density at higher energies. Many issues in cosmology like in- 
flation, dark energy, cosmological constant were investigated in the braneworld cosmological 
scenario and encouraging results were obtained. 

Although the present universe appears homogeneous and isotropic in its overall structure, 
as indicated from recent WMAP data that cosmic microwave background is isotropic up to 
1 part in 10 5 . But there are reasons to believe that it has not been so in all its evolution and 
that inhomogeneities and anisotropies played an important role in the early universe. There 
exists a large number of anisotropic cosmological models, which are also being studied in 
cosmology, due to various reasons [4J. 

Hence, it is natural to ask how these anisotropies play a role in the context of the 
braneworld scenario. The anisotropic braneworld with a scalar field is studied in [5j and it is 
shown that a large initial anisotropy introduces more damping in the scalar field equation of 
motion resulting in greater inflation. Cosmological solutions for the Bianchi-I and Bianchi-V 
in the case of the RS branemodel were studied in 6|. It is shown that for matter on the brane 
obeying the bariotropic equation of state, the anisotropic Bianchi-I, V braneworlds always 
isotropise although there could be intermediate stages in which anisotropy grows. The shear 
dynamics in the Bianchi-I brane model were studied in 7| and shown that for 1 < 7 < 2 
shear has maximum value during the phase transition from nonstandard to standard cosmol- 
ogy. The cosmological solution of field equations for an anisotropic brane in the generalised 
RS model were obtained in js| and the solution admits an inflationary era. Dissipation of 
anisotropy on the brane is explicitly demonstrated by the particle production mechanism in 
9] by considering the Bianchi-I braneworld model. Apart from this there have been studies 
lOj which show that anisotropy plays an important role in braneworld models. All the 
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previous studies have mainly focussed on anisotropy on the Randall-Sundrum braneworld. 

In this paper we consider the DGP braneworld model with anisotropy and find a solution 
of field equations. The DGP cosmological model possesses two classes of solutions; one which 
is close to standard FRW cosmology and the other is a fully five dimensional regime or a 



self-inflationary solution which produces accelerated expansion. It is also noted in [llj that 
the self-accelerating universe exhibit ghost and tachyonic like excitations. We are mainly 
interested in effects of initial anisotropy and show that it is possible to get self-inflationary 
solution in the presence of anisotropy. We also explore an acceleration condition for the 
DGP cosmological model with a scalar field dominated universe, in the presence of shear. 

II. FIELD EQUATIONS IN THE DGP MODEL 

We consider the DGP model, where our universe is a 3-brane embedded in 5D bulk with 
an infinite size extra dimension and there is an induced 4D Ricci scalar on the brane, due 
to radiative correction to the graviton propagator on the brane. In this model there exists 
a length scale below which the potential has usual Newtonian form and above which the 
gravity becomes five dimensional. The crossover scale between the four-dimensional and 
five-dimensional gravity is given by, 

where [i 2 = 87r(?4 and k 2 = 8nG5 are the constants related to 4D and 5D Newton's constants 
respectively. The Einstein equation for a five dimensional bulk is given by, 

G AB = -Ag AB + k 2 5 (T AB + S AB 8{y)), (2) 

where A is the bulk cosmological constant, Tab is the five- dimensional energy momentum 
tensor (A,B = 0,..4 ) and Sab is the 4D energy-momentum tensor and is given by, 

Sab = tab — &Qab — H 2 Gab, (3) 

where tab is the energy-momentum tensor of matter fields on the brane, a is the brane 
tension, and the last term in the above equation represents contribution from the induced 
curvature term. Assuming Z2 symmetry for the brane and using the Israel junction condition, 
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the effective Einstein equation on the brane is obtained as 12|, ll3] 



ak 4 k 4 k 4 

+ ~1 F ^ + k t U ^ + ~\ L ^ - E v ( 4 ) 

n ^ = ~\vu + 7^ T1 > + \%vT a pr aP - ^q^r 2 , (5) 
Ffiv — —-^Gp.pG P v + r^GGfa, + -q^GapG" 13 — —q^C 2 , (6) 

V = \(G,X + r w GS) - i(rC^ + Gv) - ^(G^t** - ~GV), (7) 

where is the bulk energy- momentum tensor. It is noticed that in Eq. (j4j), apart from the 
usual quadratic matter field corrections to energy momentum, there are corrections coming 
from the induced curvature term through and L^. E^ u is the projection of the bulk 
Weyl tensor on the brane representing the nonlocal effects from the free gravitational field. 
If we define the four velocity comoving with matter as [l^] , the non-local term takes the 
form, 

Efu, = -Ar 2 c [U (u^u v + i V) + P ^ + Q» u » + ■ ( 8 ) 
where = g^ u + u^u u and 

U = -±E^u" t (9) 

c 

is an effective nonlocal energy density on the brane, 

= ~ [K a K? - \h al3 V] E aP (10) 

c 

is the effective nonlocal anisotropic stress, and 

Q, = ^,X/- (11) 

c 

represents the effective nonlocal energy flux on the brane. 

At this stage we set the bulk cosmo logical constant, brane tension to zero and consider 
empty bulk, then Eq. (j3J) becomes, 

k 4 k 4 
G>„ = ~~j~^F^ u + k^U.^ + —L^ — E^ v . (12) 

We combine all the nonlocal and local bulk corrections and these can be written in an 
compact form as, 

n __ i , 



G^y — k 5 T^ , (13) 



where, 

tot 1 11 

= Ta^v-v + + —^L^y — jjEfiv ■ (14) 

/I H K 5 

Following [15J the effective energy density, pressure, anisotropic stress, and energy density, 
for the DGP case can be calculated and are obtained as, 

1/2 1 \ 4r 2 

Ptot = J^\P 2 - J^ G ooP + —,G 2 00 j + -j£U, (15) 

If 2 2 1 \ 4r 2 

Ptot = ^[P(P + 2p) - — C o(p + p) - — Giip + —G 00 (G 00 + 2G U ) + ^tjU, (16) 



4r 2 
4r 2 



< = -rf^^ (17) 



C = jfQ»- ( 18 ) 

It is assumed that the total energy-momentum tensor and the brane energy-momentum 



tensor are conserved independently [151 ] . It can be seen from Eqs. (IHl fTTl) that results 
depend on the crossover scale r c of the theory, which is a typical feature of the DGP model. 
When r c goes to zero the effective nonlocal energy density, anisotropic stress, and energy 
flux diverge. 

III. DGP COSMOLOGY ON A BIANCHI-I BRANE 

We consider an anisotropic Bianchi-I brane model with the metric given by, 

ds 2 = -dt 2 + a 2 (t)(dx i ) 2 , (19) 

and is covariantly characterized by jf| 

D„f = 0, A li = = u /t , Q M = 0, R; u = 0, (20) 

where is the projected covariant spatial derivative, / is any physical defined scalar, 
is the four acceleration, is the vorticity and R* is the Ricci tensor of the three surface 
orthogonal to u^. 

The conservation equations takes the following form: 

p + 6(p + p) = (21) 

D v P ia> = (22) 



U + IQU + a^P^ = (23) 

where the dot denotes the u u V v and G represents the volume expansion rate and a^ v is the 
shear rate. The Hubble parameter for the Bianchi-I metric is given by Hi = ^ and one can 
define the mean expansion factor as S = (aio^a^) 1 / 3 , thus 

e = 3H = 3^=Y i H i . (24) 

i 

The Raychaudhuri equation for the Bianchi-I metric on the brane is obtained as, 



1 -k 2 

— W + cx p a au = 

3 2 



1 /„ /„, K\ 2 \ 2 1 



3 # +7^ +-p(p + 2p) 



12^ V V #7 J 4 



1 A. /„o X 



3 ^ + — (p+^ + ^p 



2^2 11 £2 
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(25) 



and Gauss-Codacci equations are 

<V + 0(7^ = 4r c 2 P M , (26) 

§e> + ^ = $(s(* + *)-^)%4#r. (27) 

It is noticed that there is no evolution equation for P^, indicating the fact that in general 

the equations do not close on the brane and complete bulk equations are needed to determine 

the brane dynamics. There are bulk degrees of freedom whose impact cannot be predicted 

by brane observers. Hence, the presence of P^ v in Eq. (|26|) does not allow us to simply 

integrate to get the value of shear as in general relativity. We can overcome this problem 

by considering a special case where nonlocal energy density vanishes or is negligible 5J i.e. 

U = 0, which is often assumed for FRW branes and leads to conformally flat bulk geometry. 

Under this assumption conservation Eq. (123]) leads to o^P^ = 0; this consistency condition 

implies a condition on evolution of P^ v as a^ v P^ v = Ar^P^P^, folowing from Eq. (1261) . As 

there is no evolution equation for P^ v on the brane ^, [ijj], this is consistent on the brane. 

Also, one should check that the brane metric with U = leads to a physical 5D bulk metric. 

This has to be done numerically as the bulk metric for the Bianchi brane is not known and 

is beyond the scope of the present paper. 

Now Eq. (1261) can be integrated after contracting with shear {5]] and gives 

6Y 2 

a^a^ = 2a 2 = — , E = 0. (28) 
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FIG. 1: Evolution of scalar factor with and without anisotropy, for r c = 2 and E = 0.001 for 
e = +1. The dashed line shows the behaviour in the isotropic DGP model and the thick line 
corresponds to Bianchi-I DGP brane. It is observed that the universe enters the self-accelerating 
phase much later compared to the isotropic case. 

A. Friedmann equation for Bianchi-I Brane 

The generalised Friedmann equation in DGP cosmology for the Bianchi-I brane is ob- 
tained using metric (1T91) and Eqs. ( 1271 [281) as 
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H + ~ & ~ ~Wl {H2 + - & = ~ p - (29) 



161 ] for isotropic universe 



The above equation reduces to the DGP Friedmann equation of 

when S = 0. Here, e = ±1 corresponds to two possible embedding of the brane in the bulk 

space-time. We can obtain the Bianchi-I equation for general relativity from Eq. (129]) under 

the condition 

K Y 2 ?n 2 
(tf 2 + ^)-^»|p (30) 



which matches the crossover scale set by [3J, but in this case it also depends on the shear. 
Also see that (1301) reduces to the condition of Deffayet for recovery of standard cosmology, 
in the absence of shear. 

It can be seen from the Eq. (1291) that we recover the Friedmann equation of jlj] when \i 
goes to infinity. It also corresponds to the fully five dimensional regime, as H goes linear to 
p. This can be written equivalently as, 



K V, 2 9n 2 

(tf 2 + f;) -§««=§- OH 



which is the condition to recover the fully 5D regime from (1291) . 
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FIG. 2: Evolution of a scalar factor with different values of E, thick line (£=1), line (£=0.1) , 
dashed line (£=0.001), for e = +1 with r c = 2. The figure shows that the higher value of £ implies 
entering in to the self-inflationary phase much earlier than the smaller value of £. 

B. Solution in the late universe 



One of the features of DGP cosmology is that either it enters a fully 5D regime or a 
self-accelerating phase, in the late universe, depending on the value of e. Equation (1291) can 
be written as (with K = 0) 

2 




5! 

S 6 ' 



(32) 



The above equation can be expanded under the condition // 2 /? r c ^ 1 [following from (13TT)]. 
for the e = 1 case at zeroth order we get, 



Tj2 I 

" rl S 6 



(33) 



which can be compared with the corresponding equation of |16l . Il7j and in the absence of 
E matches their result. Also, for S — > oo as expected the shear term damps and Eq.( l33l) 
conforms with 



161 ]. Then the solution of Eq. (1331) is obtained as 



S 3 = r c S Sinh(-) 



(34) 



which shows self-acceleration in the late universe; this feature is already shown by Friedmann 
Eq. (!33l . Hence, it indicates that the DGP model can lead to self- inflation even in the 
presence of anisotropy. Figure 1 shows behavior of the scale factor for the anisotropic and 
the isotropic case. Figure 2 and 4 shows growth of S for different values of £ and r c . 



S 



Similarly Eq. ([3*2]) can be expanded for the e = — 1 case and it gives 

H 2 = |J, (35) 
which shows that the expansion rate is dominated by shear term. The solution is, 

S 3 = Sl + 3Z(t-t ). 

The expansion rate in the present case is slower in comparison to standard cosmology and 
the result matches results of [s| for the RS type brane model . This is expected in the case 
of e = — 1 which corresponds to the fully 5D regime and is the same as the RS model when 
quadratic term dominates. Figure 3 shows behavior of S for different values of S. To discuss 
the early universe scenario, we rearrange Eq. ( 129]) as follows: 

2 

V2 

3 \2r,. V/'V V ' <■ •>■■■' ■> ^' ^ 

At high energies the above equation can be expanded in terms of fJ?pr 2 3> 1. Thus at zeroth 
order we get the equation of the Bianchi-I model of the general relativity. 

C. Acceleration conditions 

Next we get the acceleration conditions in an anisotropic DGP brane model. Consider 
Eq. (1321) up to first order for the e = +1 case. We obtain the Friedmann equation, 

The accelerated expansion in standard cosmology is given by H + H 2 > 0. Thus using Eqs. 
( [""!]) and eqn fl371) the acceleration condition for the present case becomes 




and this gives, 



H + H 2 ~ -/i 2 (p + p) + i + ^ 2 p - 2^ > 0, (38) 
P 1 2E 2 

Similarly for the e = — 1 case, expanding (132]) and considering the next order, the Friedmann 
equation is, 

..4 2 y2 

H 2 = + |f • (4°) 
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FIG. 3: Evolution of scalar factor for the case e = —1, with different values of E. Dashed line ( 
E=l), thick line (£=0.1), line (£= 0.001). 

Notice that the relation between the Hubble radius and energy density is linear, a feature 
of brane 
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181 ]. which is referred to as the fully 5D regime. The acceleration condition, 



H + ~ -^p(p + p) + - 2^ > 0, (41) 



implies 



2 P 6S2 

p< -Y~W^p- (42) 

Therefore, one can see that the condition is modified compared to the standard cosmology 



and our results reduce to 13] in the absence of anisotropy. 



Next we consider the inflationary phase driven by a scalar field, with energy density and 
pressure given by p = \<\> 2 + V(<f) and p = \§ 2 — V((f)), respectively. The Klein-Gordon 
equation for a scalar field is + 3H(j) + V'((f>) = 0. Using Eq. (!39l) the slow-roll condition 
for inflation can be obtained for e = 1 case as, 

^+2^-f- < 43 > 

It can be seen that the condition depends on the anisotropy and in the ordinary DGP 
conditon is recovered in absence of £ . Similarly for e = — 1, using eqn (H2l) the slow roll 
condition is obtained as, 

P < -*V + 2j3V> + ™ i (44, 
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FIG. 4: Evolution of scalar factor in the e = +1 case, with different values of r c = 2,1,0.1 for 
S = 0.001. The figure shows that scale factor grows faster for small values of r c , compaired to 
higher values. 

IV. CONCLUSIONS 

To summarize, we considered anisotropic effects in the DGP cosmological scenario and 
found a solution to the corresponding field equation. We obtained the Friedmann equations 
with the Bianchi-I brane and two branches of solutions in the DGP model are considered. 
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The evolution of a scale factor in the case of the late universe is studied and an acceleration 
condition for inflation is derived. 

It is observed that, for the e = +1 case, there exists a self-inflationary solution which 
leads to accelerated expansion in the late universe, even in the presence of anisotropy. The 
evolution of a scale factor is given in Fig. 1, which shows that the presence of shear does not 
stop the universe from entering a self-accelerated phase, but slows it down. The evolution 
scale factor in the anisotropic DGP model depends on the cross-overscale and the shear. It 
is noted that, when shear is unity(S = 1), the evolution of a scale factor coincides with the 
isotropic DGP model. Fig. 2 shows dependence of a scale factor on the shear and the higher 
the value of £ the faster the universe enters into a self-accelerated phase. Fig. 4 shows the 
behavior of a scale factor with typical values of a crossover scale r c and shows that the higher 
the r c the slower the acceleration. Hence, in the DGP model shear and r c behave oppositely 
to each other. The behavior S for the e = — 1 branch of the solution of the Friedmann 
equation is shown in Fig. 3 and a scale factor grows faster for higher values of anisotropy. The 
scale factor grows as S oc t 1//3 which implies slower expansion than the standard cosmology 
and is similar to expansion rate that of Randall-Sundrum case. The Friedmann equation 
in this case is similar to the RS Friedmann equation with quadratic energy density. The 
acceleration conditions and slow-roll conditions for inflation are obtained and they depend 
on the shear. 

Finally, our results coincide with those of Deffayet whenever £ is absent and reduce 
to the Bianchi-I of the general relativity under condition ( 1301) . Hence, the presence of initial 
anisotropy does not adversely affect the features of the DGP model. It is possible to have 
an anisotropic brane in the DGP model and still get a self-accelerating universe in the late 
time, as the anisotropic term appearing in fl33l) is not like energy density. Also, for a large 
value of S (S — > oo), the anisotropic term tends towards zero and we recover the isotropic 
case. 

It can be noticed that in the present work, ghost like instabilities in the self-acclerated 
branch, as pointed out by 11] . do not appear. This may be due to the fact that we are using 
special boundary conditions. So, it is interesting to study the self-accelerated branch with a 
ghost term, by solving the full brane-bulk system and considering more realistic boundary 
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conditions, which is beyond the scope of the present paper. 
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